function [d]=ex11_2_3(t,y)
    d=[y(2);(1-y(1).^2)*y(2)-y(1)];
end

% [T,Y]=ode45(@ex11_2_3,[0,20],[2,0]);